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Neutralino relic density from direct non-equilibrium production 
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We review the calculation of the LSP relic density in alternative cosmological scenarios where contributions 
from direct decay production are important. We study supersymmetric models with intermediate unification 
scale. We find concrete scenarios where the reheating temperature is of order one GeV (Mi ~ fO 12 ) and below. 
If the case that this reheating temperature is associated to the decay of oscillating moduli fields appearing in 
string theories, we show that the LSP relic density considerably increases with respect to the standard radiation- 
dominated case by the effect of the direct non-thermal production by the modulus field. The LSP can become a 



good dark matter candidate (0.01 - 0.1 > Qh 2 < 0.3 - 1) for Mi 



10 1 



10 14 Gev and ~ 1 - 10 TeV. 



1. Introduction 

In a series of works various implica- 

tions for an intermediate unification scale in su- 
persymmetric models have been considered. In 
particular, it was pointed out that the lightest 
neutralinos (LSP)-nucleon cross sections are very 
sensitive to the variation of this scale. For in- 
stance, for Mi ~ 10 12 - 10 14 GeV the cross 
sections are enhanced and a large region of the 
parameter space of the minimal supersymmet- 
ric standard model (MSSM) is compatible with 
the sensitivity of the current and near future di- 
rect detection dark matter experiments (DAMA, 
CDMS, EDELWEISS). In fact, such results can 
also be obtained within the GUT unification sce- 
nario, i.e., Mi ~ 10 16 GeV, but it requires tuning 
the available parameter space: either large tan/3 
(tan/3 > 25) or a specific choice for non-universal 
soft SUSY breaking terms so that the value of /i 
term is reduced and hence a larger Higgsino com- 
ponents of the LSP is obtained, which is essential 
for enhancing the LSP-nucleon cross sections ||. 

One problem is than in GUT unification sce- 
narios with standard Big Bang (BB) cosmological 



expansion where the LSP is produced in a ther- 
mal environment dominated by radiation, such 
large cross section leads to a very small LSP relic 
density (outside the cosmological and astrophys- 
ical interesting region 0.1 < Vlh 2 < 0.3). In this 
scheme there is a generic anti-correlation between 
the LSP relic density and its quark or nucleon 
cross sections: one usually obtains 



n x h 2 



C 



(i) 



"The work presented here was made in collaboration with 
S. Khalil and C. Munoz, see. Ref.jl] 



where a ann is the annihilation cross section for 
neutralino pair, v is the pair relative velocity, and 
(..) denotes thermal averaging. Scattering (cr x _ g ) 
and annihilation cr ann LSP cross sections can be 
related by a cross symmetry. Therefore, a large 
scattering cross section u x - q leads generically to 
a large annihilation cross section and a small relic 
density as Eq. 1 suggest. Indeed, in these models, 
it was observed that LSP-proton cross sections of 
order 10 -6 Pb correspond to too low relic densi- 
ties (<~ 0.005- 0.01 ). 

We want to address here the relic density ques- 
tion: in the context of diverse non-standard cos- 
mological scenarios, what is the effect of chang- 
ing the unification scale from GUT scale (~ 10 16 
GeV) down to intermediate scale (~ 10 12 GeV) 
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and whether we still have any anti-correlation be- 
tween il x h 2 and a x m . 

2. The standard BB scenario 

The LSP relic density is given by 

$l x h — Pxl Pcriti 

where p x is the LSP mass density and p cr u = 
3H$/8ttG ~ 1.9 x 10~ 29 h 2 g/cm 3 . The evolution 
of the number density n x = p x /m x can be de- 
scribed to a precision sufficient for our purposes 
by the Boltzmann equation 



dn x 
~~dt 



3Hn x = 



l v) [(r 



(2) 



where H is the Hubble expansion rate. In the 
standard calculation, it is assumed that in the 
early universe (T ^> m x ) the LSP density was 
initially equal to its equilibrium value n x = n x q . 
As the universe expands the LSP density traces 
its equilibrium value. When the particle becomes 
non relativistic ( T < m x ) the number density of 
x's drops exponentionally. The annihilation rate 
r x = (<r x nn v)n x becomes smaller than the expan- 
sion rate, T x < H. At this point, the LSP decou- 
ples and fall out from being in equilibrium with 
the radiation, their number remains constant af- 
ter that and its density is only diluted by the 
universe expansion. 

As long as the LSP is non-relativistic, we can 
expand the averaged cross section (a x nn v) as fol- 
lows 



(a a x nn v) = a s +a p (v 2 ), 



(3) 



where a s is the s-wave contribution at zero rel- 
ative velocity and a p contains contribution from 
both s- and p- waves. The relic density solution 
to the Boltzamn equation is accurately given by 
the approximate formula (c.f. Eq.l) 



8.8 x lQ- n GeV~ 2 
y/g*(a s /x F + 3a p /x 2 F ) 



(4) 



where the freezing-out temperature x F = m x /Tp 
can be estimated as 



xp ~ In 



c(a s + 6(Xp/xp) 
\/g*x F 



(5) 



In these formulas <?* is the effective number of 
relativistic degrees of freedom at T F and c is con- 
stant of order 1CT 1 jjl). 

3. Alternative scenarios 

In alternative cosmological scenarios, it is often 
assumed the existence of one (o more) epoch be- 
fore the radiation-dominated era of the universe 
where the energy is dominated by coherent oscil- 
lating fields, the so called reheating era. A reheat- 
ing process is associated by example with the final 
stage of Inflation. It can be also, independently or 
not, associated with the oscillations of moduli or 
Polony fields appearing in string theory after flat 
directions acquire a mass from SUSY breaking. 

The details and physical origin of these reheat- 
ing are largely unimportant. Its effects can be 
characterized by a reheating temperature Trh 
which is conventionally defined (see the definition 
below) as the temperature at which the oscillat- 
ing field energy ceases to dominate the cosmolog- 
ical evolution and starts the radiation dominated 
epoch. 

The standard BB scenario presented in the pre- 
vious section can be seen as an special case when 
Trm is too much higher than Tp (of order 10 9-10 
GeV ). In this case, the reheating epoch has no 
relevance in the final output of the relic density. 

A problem could appear if the decay of the 
oscillating field occurs after or during the pri- 
mordial BBN nucleosynthesis. This is drastically 
solved imposing from the onset that the reheating 
temperature is larger than ~ 1 MeV. 

Finally, in interesting intermediate cases, a rel- 
atively low value of T^h such that Tp > Tr/j can 
have qualitative and quantitative implications on 
the predictions of the relic abundance of the LSP 
as discussed in Ref. The main reason for 

this is that now the decoupling of the LSP par- 
ticle occurs in an expanding universe dominated 
by matter instead than radiation and that has 
effects. We will see below that theories with in- 
termediate unification scale predict low reheating 
temperature ( < GeV and below) and they can 
be considered as interesting scenarios. 

It is useful to use the following expression which 
can be interpreted as a quantitative definition of 
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the temperature Trh and which connects it to 
the decay width of the oscillating field: 



^9,(T RH ) 
45 



-1/2 



T 2 
_ RH 

M P ' 



(6) 



Given any concrete theory, the reheating temper- 
ature can be related to the their physical param- 
eters inverting Eq.^| for Trh as we can see on the 
following. 

Let us consider first an inflation scenario. The 
inflaton ((f)) decay width which can be writ- 
ten, in some some supersymmetric hybrid infla- 
tion which has been considered |17l, as 



U = 



1 M 



(0) 



where is the inflaton mass. Here Mf is the 
mass of the particle / that the inflaton decay to 
(i.e. a right handed neutrino or a sneutrino). Mf 
should be less than the inflaton mass to allow 
for the decay (j> — ► //. Finally (</>) is the vev 
of the inflaton field which is of the order of the 
unification scale. 

In Ref. |Q it has been shown that a intermedi- 
ate unification scale of order Mj ~ 10 n ~ 12 GeV 
is favored by inflation. In this case, recalling that 
the inflaton mass is constrained by 

m < Mf/Mp, 

one obtains ((0) ~ 1.7x 10 16 GeV |l7)) an inflaton 
mass of the order ~ 10 2 GeV. For the 'stan- 
dard GUT scenario' with Mi ~ 10 16 one obtains 
however ~ 10 13 GeV . From Eq.^) we find 
that Trh — C(l) GeV; a value much smaller than 
the values quoted in the standard GUT scenario 
(Trh ~ 10 11 ). This reheating temperature is of 
the same order or lower than the typical freeze- 
out temperature Tp ~ m x /20 ~ 0(1 — 5) GeV. 

A detailed analysis of the relic density with a 
low reheating temperature has been considered in 
Ref. JTj| with different possible generic scenarios. 
The first scenario corresponds to a case where the 
LSP is never in chemical equilibrium either before 
or after reheating. In the second the LSP's reach 
chemical equilibrium but they freeze out before 
the completion of the reheat process Tp > Trh ■ 



Each of these scenarios leads to different qual- 
itative and quantitative predictions for the relic 
density which can become quite different from the 
standard computation we summarized in Eq.(|]). 

CASE A) In case of non-equilibrium produc- 
tion and freeze out at the early times and sup- 
posing that the number density of the LSP n x is 
much smaller that n x q , the Boltzmann equation 
(|^) can be approximated and solved. In this case 
one gets |L3] 



Qh 2 



2.1 x 10 4 A, (a. + 2E 



(7) 



where 



A, 



'.£_ 

4 



g*(T 



RH) 



10 



.3/2 



10 



g* (T* 



(10 3 T 



RH) 



(8) 



and where g is the number of degrees of free- 
dom of the LSP and T* is the temperature at 
which most of the LSP production takes place, 
it is given by T* ~ 4m x /15. As we can see this 
flh 2 is proportional to the annihilation cross sec- 
tion, instead of being inversely proportional, as in 
Eq.(|4|). However, the assumption that n x « n x q 
leads to a sever constraint on the annihilation 
cross section IrL 



a s < a s and a p < a p where 



io- 



10 J and 



10" a - 10" 



So 



that for large neutralino cross sections, which we 
are interested in, eq.(0) can not be applied. 

CASE B) With large annihilation cross section 
the LSP reaches equilibrium before reheating as 
discussed in Ref.|l3| and its relic density Vth 2 is 
very close to the one obtained by using the stan- 
dard calculation of relic density in Eq.(|^). In this 
case, Tp is obtained by solving a different equa- 
tion as before 11131 



xf = In 



/ 5/2,5 3/2 

G I a s Xp + -a p Xp 



(9) 



where C is an unimportant constant now. The 
relic density flh 2 is given by 



, _ gl ,2 (T RH ) 2.3 x 10- n T 3 H GeV- 2 



g*(T F ) m x (a s x F A + 4a p x F 5 /5) 



(10) 



In this case the relic density is again inversely 
proportional to the annihilation cross section as 
in Eq.(0). Moreover, it has a further suppression 
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due to the very low reheating temperature effect. 
The typical value for the relic density in this case 
for a ann ~ 1(T 6 GeV~ 2 is Vth 2 < 1CT 3 . So it is 
even worse result than the BB standard compu- 
tation's one. 

4. Our Model 

In superstring theory, which is the first moti- 
vation for the idea of considering an intermediate 
scale ||, there are many moduli fields which ac- 
quire masses from SUSY breaking effects. Their 
masses are expected to be of order the gravitino 
mass but their couplings with the MSSM matter 
content are suppressed by a high energy scale (l5) . 
The moduli ((f)) decay width will now be param- 
eterized now as 



where Mj is the unification scale (~ 10 12 GeV) 
which acts as an effective suppression scale. 

The reheating temperature is obtained using 
eqs. (§) and (0) 

T RH _ 10 14 GeV / \3/2 / 10.75 V 

5 MeV ~ Mi V100 GeV J \g*(T RH ) J 

where, for example, g* — 10.75 for T ~ 0(1 — 10) 
MeV. This reheating temperature is shown as a 
function of the modulus mass in Fig. [I] for dif- 
ferent values of Mj. The request that the mod- 
ulus mass is larger than ~ 100 — 500 GeV in 
order to allow for kinematical decays into neu- 
tralinos of suitable mass m^o ~ 50 — 200 GeV, 
limits in practice the reheating temperature to 
be above ~ 3 GeV for the lowest scale on con- 
sideration Mi = 10 11 GeV. This has important 
consequences for the relic density computation. 
As discussed in the introduction of this section, 
we need a temperature smaller than the typical 
Tp ~ m^o/20 ^3—10 GeV, in order to in- 
crease the relic neutralino density. This implies 
that the scale Mi ~ 10 11 GeV is in the border 
of validity. On the other hand, for larger val- 
ues we can obtain very easily interesting reheat- 
ing temperatures. For example, for Mj = 10 12 
GeV we have T RH > 0.3 GeV. For the highest 
scale with an interesting phcnomenological value 



of the ncutralino-nucleus cross section, in the case 
of universality and moderate tan/3 j|, Mi = 10 14 
GeV, the lowest value of the reheating tempera- 
ture corresponds to Trh ~ 6 MeV. Larger val- 
ues of the scale, Mi > 10 15 GeV, producing also 
a large cross section, are possible in D-brane 
scenarios since non-universality in soft terms is 
generically present |I]J^]. In this case the con- 
straint Trh > 1 MeV from nucleosynthesis can 
be translated into a constraint on the modulus 
mass > 140 GeV. 




Log 10 m t (GeV) 



Fjigure 1. The reheating temperature Trh as 
a function of the modulus mass for different 
values of the effective scale Mi = 10 11 , 10 14 , 10 16 
GeV. The horizontal at O(l) MeV is the nucle- 
osynthesis limit. 



When considering the direct decay of the mod- 
uli fields into LSP (and therefore, an additional 
quantity of non thermal LSP production) the evo- 
lution of the cosmological LSP abundance be- 
comes more complicated. We have to consider 
now the coupled Boltzmann equations for the 
LSP, the moduli field and the radiation [p^|JT3]| . 
The detailed equations used by us appear explic- 
itly in Ref.(0). These equations depend on m x , 
the LSP mass, Nlsp, the averaged number of 
LSP particles produced in the decay of the mod- 
ulus field. 

Let us discuss qualitatively the solution follow- 
ing the arguments used in ref. [^5|. For a Trh 
higher than Tp the relic density will roughly re- 
produce the usual result given by eq. (H). How- 
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ever, for the interesting case for us when Tun is 
lower than Tp, neutralinos produced from modu- 
lus decay are never in chemical equilibrium, un- 
like the thermal production case reviewed in Sec- 
tion 2. As a consequence, its number density al- 
ways decreases through pair annihilation. When 
the annihilation rate (cr??"u)n^o drops below the 
expansion rate of the universe, H, the neutralino 
freezes out. Then the relic density can be esti- 
mated as [|l|,|l5[ 



n^oh 2 

Al 



3 m? 



2 (2^2/45) g* Tl H (a 



Pel 



43.) 



This result is valid when there is a large number 
of neutralinos produced by the modulus decay. 
When the number is insufficient, they do not an- 
nihilate and therefore all the neutralinos survive. 
The result in this case is given by 



n^oh 2 = 

Al 



3 N^o m x T 2 M 2 



(27r 2 /45) g* T^ H m<j, p c /s 



(14) 



The actual relic density is estimated [|15| as the 
minimum of dl3| ) and (|l4|). 

Now we can apply the above equations to our 
case with intermediate scales. Using eqs. (O) and 
12|) , we can write expressions (Of) and (F[j) as 



il^oh 2 cx 



Mim^o 
ai 



[a^vW^gl^ 



(15) 



constant relic neutralino density fl^o h 2 as a func- 
tion of and N~ x o are shown, for fixed values 
of Mi and m^o. In particular, we consider the 
cases Mi = 10 12 , 10 13 , 10 14 GeV, with m x o = 100 
GeV. The corresponding reheating temperatures 
can be obtained from Fig. 1. Note that whereas 
many values of iV^o correspond to a satisfactory 
relic density for Mj = 10 12 - 10 13 GeV, for the 
case Mi = 10 14 GeV only a small range works. 
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(16) 



From these equations we can see that even with a 
large annihilation cross section, (a^} n v) ~ 10~ 8 

GeV -2 , we are able to obtain the cosmologically 
interesting value f^o/i 2 



'XV 

Mr ~ 10 13 GeV we obtain it when N?o 



1. For example for 

using eq. (|15|). 

Previous expressions are useful but approxi- 
mate. We have solved numerically the Boltzmann 
equations in several cases. In Fig. || we show 
in detail the results of the numerical computa- 
tion. We have used a particular expression for 
the annihilation cross section introduced which is 
explained below. In this figure, the contours of 



Figure 2. 
as a function of and N for m x 



Relic Neutralino density Qti 2 contours 
= 100 GeV and 

values of the intermediate scale Mi = 10 11 ,10 14 
GeV. 



Let us finally finish this section with a few 
words about the neutralino annihilation cross sec- 
tion. As it is well known, in most of the parameter 
space of the MSSM the neutralino is mainly pure 
bino, and as a consequence it will mainly annihi- 
late into lepton pairs through t-channel exchange 
of right-handed sleptons. The p-wave dominant 
cross section is given by ]l2] , ^3| 



(a~o v) 



1 



1 



iira ' 



(17) 
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where xr = m~ /m 2 and a' is the coupling con- 

Ir in 1 x{ 1 b 

stant for the U(l)y interaction. Taking mj ~ 



m^o ~ 100 GeV, (cr?o n w) in eq. ( |17| ) becomes of 

the order of 10~ 9 GcV~ 2 or smaller. Using eq. (Q) 
an interesting relic abundance, O^o/i > 0.1, is 
obtained. 

However, in Refs.^jQ], we obtained that for 
Mi ~ 10 12 - 10 14 GeV the neutralino LSP ac- 
quires a sizeable Higgsino component, as we com- 
mented before this is essential for enhancing the 
LSP-nucleon cross section. For a Higgsino-like 
LSP the dominant annihilation channel is the de- 
cay into a W boson pair. An upper bound on the 
cross section can be obtained 



<= 



1 (l-x w ) 



3/2 



2 m 2 {2-xwf 



(18) 



where ai is the SU (2) constant and xw = 
m^/m 2 ,. It is this last expression for the cross 
section the one used in the numerical example 
presented in Fig.[j]. For a LSP with m x ~ 100 
GeV, the formula gives a maximum ~ 2 x 10 -8 
GeV -2 . In the case that the LSP have sizeable 
components of both gaugino and Higgsino one ex- 
pects smaller cross sections. 

5. Conclusions and Comments 

In conclusion we have analyzed the LSP relic 
density in supersymmetric models with interme- 
diate unification scale, we find that the reheat- 
ing temperature is of order 1 GeV (Mj ~ 10 12 ) 
and below down to the nucleosynthesis thresh- 
old. We have shown that in generic low Trh sce- 
narios without non-thermal direct production is 
expected that the LSP relic density is too low. 
As an alternative, we have presented a concrete 
scenario where the reheating temperature is asso- 
ciated to the decay of oscillating moduli fields ap- 
pearing in string theories, we show that the LSP 
relic density considerably increases with respect 
to the standard radiation-dominated case by the 
effect of the direct non-thermal production by the 
modulus field. The LSP can become a good dark 
matter candidate (0.01-0.1 > Vth 2 < 0.3-1) for 
Mi ~ 10 12 - 10 14 GeV and ~ 1 - 10 TeV. 

Let us finally remark that the numerical value 
of N^o is in general model dependent. This 



was discussed in the context of supergravity in 
ref. In this particular case both values 

N^o ~ 1 and N^o ~ 10~ 3 — 10 -4 are plausible, de- 
pending on the characteristics of the supergravity 
theory under consideration. 
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